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1 Introduction 



One of the most important recent developments in the field of heavy-quark physics has 
been the formulation of the static limit for a heavy quark as an effective theory 
The so-called Heavy Quark Effective Theory (HQET) has put the description of heavy- 
hadron physics on a QCD-related and model-independent basis. Aside from this theo- 
retical progress there have been numerous phenomenological applications of this idea, in 
which the model dependence has completely disappeared or is substantially decreased 0]. 

Most of these applications deal with systems involving a single heavy quark. In HQET 
both the particle and the antiparticle numbers are separately conserved, and the appli- 
cations considered so far deal mainly with the one-(anti)particle sector of HQET. Thus 
the one-particle sector of HQET is already explored in some detail, including leading and 
next-to-leading QCD corrections. 

However, there have also been first attempts to deal with the two- (anti) particle sector 
and the particle-antiparticle sector |^. Here it turns out that complex anomalous 
dimensions are obtained by a naive analysis. In HQET the anomalous dimensions in 
general depend on the velocities of the heavy quarks and the imaginary parts of the 
anomalous dimensions in the two-particle sector diverge as 1/ \J {vv'Y — 1 for v v', if v 
and v' are the velocities of the two heavy quarks. Subsequent investigations have shown 
that this singularity is related to the long-range part of the quark- ant iquark potential 
and that one may remove this divergence by a suitable definition of the multiparticle 
states of HQET. In a basis consisting of the redefined states one then has real anomalous 
dimensions. 

The physical origin of these phases is the same as that of the well-known Coulomb 
phases. The one-gluon (one-photon) exchange potential decreases too slowly, leading 
to the well-known singularities of the form l/\v\, where v is the relative velocity of the 
outgoing coloured (charged) particles. From its physical interpretation as well as from the 
fact that the divergent phase is a property of the states, it is clear that these divergent 
contributions are a long-distance effect belonging to the infrared dynamics of the state 
with two heavy quarks. 

Heavy quarkonium states have to show up as bound states in the particle-antiparticle 
sector of HQET. In such a state, however, the two velocities of the heavy quarks differ by 
an amount of only Aqcn/fn (^ is the heavy-quark mass) and hence we would rather like 
to switch to a description in which the two heavy-quark velocities become equal. This 
limit, however, cannot be taken in a naive way, since the static limit of HQET does not 
reproduce the above phases, which are related to the potential between the two heavy 
quarks and hence to the binding mechanism of the quarkonium state. 
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It turns out that this hmit may be performed if the evolution of the states is determined 
not by the static Lagrangian, but rather by a Lagrangian consisting of the static HQET 
part and the first subleading spin symmetric contribution to the Lagrangian, i.e. the 
kinetic energy operator. In other words, the presence of the divergent phase forces us to 
go beyond the static limit, if we want to consider heavy- quarkonium states, in which the 
velocities of the heavy quarks differ only by an amount of order Kqcd/^- This in turn 
introduces a mass dependence into the lowest-order dynamics such that the inverse "Bohr 
radius" A, which sets a small non-perturbative scale, is in general not independent of m. 
For this reason we lose heavy- flavour symmetry, but spin symmetry is still present, since 
the kinetic energy operator is still spin symmetric. 

In the present paper we set up an effective theory approach to processes involving 
heavy quarkonia, which is based on the heavy-mass limit of QCD. We focus on inclusive 
annihilation- type decays for which we perform a systematic expansion in powers of A/m, 
up to logarithms, which may be accessed via the renormalization group. 

The approach proposed here has many features common to the one of Bodwin, Braaten 
and Lepage (BBL) P]. Their method is based on non-relativistic QCD (NRQCD), where 
the lowest-order dynamics is determined by the Lagrangian of the Schrodinger equation 
for the heavy quarks, which is basically that obtained by including in addition to the static 
HQET part also the kinetic energy operator. Based on this they perform an expansion 
in terms of f/c, where v is the typical relative velocity of the two heavy quarks. We 
shall compare the two approaches between them as well as with previous approaches to 
inclusive quarkonium decays (for reviews see e.g. P, 0) as we go along. 

Our general strategy is as follows. We shall first discuss (section 2) the heavy-mass 
limit for quarkonia-like systems and show that the two heavy-particle velocities may be 
chosen to become equal, once the subleading kinetic energy operator is added to the 
lowest-order dynamics. This however implies that there will be no mass-independent 
static limit for the quarkonia, in other words, the "binding energy" A will depend on the 
heavy mass in a non-trivial way, but it is still a small scale A <C m, so that an expansion 
in powers of A/m is useful. 

The first step to access the inclusive annihilation decays is to separate long and short 
distances, where the distance scale is set by the Compton wavelength of the heavy quark. 
Technically speaking, the first step is to set up an operator product expansion for the 
inclusive decay rate of a heavy quarkonium into light degrees of freedom (section 3). In 
section 4 we write down the A/m expansion for the rate of inclusive heavy quarkonia 
annihilation up to and including A^/m^; up to this order the inclusive decay rate involves 
matrix elements of dimension-six and dimension-eight operators. 

The coefficients of this expansion may be calculated in perturbation theory, and some 
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simple examples are considered in section 5. These perturbatively calculated coefficients 
are the ones at the matching scale, i.e. at the scale of the heavy quark mass. In general 
the coefficients are scale-dependent, since the operators need to be renormalized. In 
section 6 we consider the renormalization group flow of the coefficients. Given their 
values at the matching scale one may then use the renormalization group to scale down 
to some small scale /i, at which one then tries to estimate the matrix elements of the 
corresponding operators. These matrix elements are non-perturbative quantities, but 
heavy quark symmetries restrict the number of independent parameters. This is studied in 
section 7. Finally, we apply this machinery to some simple examples in section 8. Section 
9 contains a further discussion of previous approaches as well as our conclusions. 

Some technical details are left for the appendices. 



2 Heavy-Mass Limit for Quarkonia 

We shall first set up the heavy-mass limit appropriate for heavy-quarkonia states. We 
shall start from the Lagrangian and the fields of QCD. We denote the heavy-quark field 
of full QCD by Q and define 

Q[^\x) = exp(2m vx)Q{x) = h^+\x) + Hi+\x) , (1) 

where f is a velocity (f^ = 1), which is later identified with the velocity of the heavy 
hadron. Extracting this phase factor from the full QCD field Q removes the dominant 
part mv of the heavy-quark momentum, since this phase redefinition corresponds to a 
splitting of the heavy-quark momentum according to p = mv + fc, where the residual 
momentum k is small, of the order of Aqcd- Furthermore, /i^"^-* {H^'^^) is the large (small) 
component field, corresponding to the projections 

/iW = P+QW, //W = P_QW with P± = 1(1 ± 1^) . (2) 

The small component field H^'^^ is related to the large scale m; integrating out H^'^^ from 
the generating functional of QCD Green's functions corresponds to the replacement 

H^+) =p_( 1 ) ^ni/i{+) (3) 

\2m + wDj ^ " ^ ' 

and this yields a non-local Lagrangian of the form ||^ 

L = hi;^\zvD)hM + h[^\l^P^ ( ] . ) , (4) 

Vim + IV D J 
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which still contains all orders in 1/m. However, the non-locality appearing in the second 
term of (Q) may be expanded into an infinite series of local terms, which come with 
increasing powers of 1/m. Hence one may in this way establish the desired heavy-mass 
expansion. The first few terms of the expansion for the Lagrangian are 



V 

2 



\2m 



+ iir-) h^^^i^P-{-ivD)ii:j)h[+^ + ■■■ , (5) 

V/m/ 

while the field is given by 

gi+) = (^1 + :^P-i^ + -^{-ivD)P^i^ + 0(l/m3)) h^+\x) . (6) 

The non-local expression @) is still equivalent to full QCD; in particular it is indepen- 
dent of the still arbitrary velocity vector v. In fact, the Lagrangian is invariant under 
an infinitesimal shift of the velocity 

V ^ V + 5v V ■ 5v = Q 

iD — s> —m6v . (7) 

This invariance is the so-called reparametrization invariance 0, which will play some role 
in what follows. 

The Lagrangian and its expansion in powers of 1/m (^) is the Lagrangian for a 
heavy quark. In the infinite mass limit, the quarks and the antiquarks are separated by an 
infinitely large mass gap, and hence in order to describe heavy quark-antiquark systems 
we have to introduce the antiquark field as a separate field, the Lagrangian of which is 
obtained from (^) or (^) by the replacement v —v. 

Heavy quarkonia should appear in the particle-antiparticle sector of HQET as bound 
states. The starting point of our considerations is the Lagrangian of a heavy-quark and 
a heavy-antiquark field. In the static limit we obtain from (Kl) 



where the static quark field h[^'> {h\^ ^) moves with velocity v {w). From (H) we obtain 
the equations of motion 



[IV 



D)h^+) = {iwD)h[;^ = . (9) 
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However, (H) or cannot be used to describe heavy quarkonia states. In order to discuss 
this we shall consider the matrix elements 



{A\Q{x)rQ{x)\0) (10) 

where A is a state containing a heavy quark and a heavy antiquark moving with velocities 
V and w respectively. In the static hmit this matrix element becomes 

{A\hi-'\x)Thl-\x)\0) (11) 

where the tilde denotes the static limit of the state. 

Matrix elements of this kind have been considered already in 0, ^, where the short- 
distance corrections have been calculated. It has been observed that for the naive defi- 
nition of the states, the anomalous dimension of the current h^\x)Th!^\x) acquires an 
imaginary part of the general structure 

lm7 = /(a.)^=4— (12) 



(13) 



as V — *■ w, where the function / is known as power series in as up to two loops | 

"31 10 



/(tts) = l^tts ( 1 + ^ 
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The real part of the anomalous dimension vanishes in the limit v ^ w to all orders, since 
the current h^^"* '-^ ^h\^^ is a generator of heavy-flavour symmetry. On the other hand, the 
anomalous dimensions of h^'^^'j^^h^^ and related by analytic continuation 

and hence the real part is identical. 

Thus the only problem is in fact the imaginary part given in (0), which diverges 
in the limit v w, and this is the major obstacle in taking the naive limit v w in 
(§). Keeping only this imaginary part, the solution of the renormalization group equation 
yields a phase factor 

exp{i(j){vw)) = exp { i / da—— ) ; /3(a(/i)) = /i— a(/i) (14) 

[vwY - 1 / , /^(«) ^/^ 

^ ' OLs{m) ) 



for matrix elements such as (|TT]), which is not well defined in the limit v ^ w. 

It has been pointed out in 0] that these phases are related to the Coulombic part 
of the one-gluon exchange and that the phase may be removed by a suitable definition 
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of the multiparticle states of HQET. After the redefinition of the states, the anomalous 
dimensions are real and well behaved in the limit u — >■ w. Furthermore, the phase in (JT^) 
is related to the "long-range" part of the one-gluon exchange potential and indicates the 
possibility of having bound states in some of the quark- (anti) quark channels. 

Consequently, the phase appearing in ([l^) is an infrared contribution which should 
be contained in the dynamics of the effective theory. In other words, if we want to 
describe heavy quarkonia, the phase factor is related to the binding mechanism, which is 
an infrared effect. In a state such as a quarkonium, the two velocities differ only by a 
small amount of order 1/m which is a hint that we need to go beyond the static limit to 
describe quarkonia states. 

In order to see how higher-order terms in the Lagrangian cure the problem, we use 
the reparametrization invariance (|^. The expressions 

W iD 

V = v ^ and W = w H 

m m 

hi-"' = (l + ^) and = (l + ^) (15) 

are invariant under the reparametrizations v + Sv and w + 6w, where this is true for 
the second line only to order 1/m. This observation has been used in []TU| to obtain, for 
heavy-heavy currents the renormalization of the subleading terms from the leading ones. 
Using reparametrization invariance, we have 



exp 



'Mvw))hi+'>Th^s2 ^ C^exp(#(VW))r/iL~2 . (16) 



If we now consider the limit w — > w we have also V W, which yields in the phase factor 
the formal expression 



exp(i0(VW))-.exp(z(/.(V^)) = exp ^-^_= J da^\, (17) 

where we have only kept the singular term in the last step. 

Hence the divergent phase factor may be shifted from the velocities into the residual 
momenta, if the equation of motion of the heavy-quark field is not the one given by (P), 
but rather by 



|(V^-1)/.W = and f (w^-l)/.L-2 = 0. 
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In other words, the Lagrangian leading to these equations of motion is the combination 



= h[^\rvDW^ + h^^'^^h^^^ - /.L-)(.^Z^)/.L-) + /^l-^^e^ , (19) 

where we have replaced w — > —w in the last step. This expression for the Lagrangian is 
reparametrization-invariant; it is also the same invariance which ensures that the static 
and the kinetic terms of the Lagrangian renormalize in the same way. Now one may per- 
form the limit f — > w without encountering a problem in the short- distance contributions 
as in (|T2p. Of course, the phase (|T3p has not disappeared: it will show up as a singularity 
in the residual relative momentum, which will be generated by the infrared dynamics as 
given in (O) . This also means that it is a long-distance effect and may be absorbed into 



the states as discussed in [g]. 

However, if it is ([T9|) that determines the leading term of our expansion, we will not be 
able to perform a strict infinite-mass limit, since now Cq depends explicitly on the heavy 
mass. In the case of heavy quarkonium, non-perturbative effects will generate binding of 
the two heavy objects that will introduce a small scale A, which now in general depends 
on the heavy mass. 

The bottomonium and the charmonium are far from being Coulombic systems and 
there is no obvious reason why in the heavy-mass limit a heavy quarkonium should become 
Coulombic; still the case of a Coulombic system is instructive. Neglecting any running 
of a, the size of such a Coulombic system is Rsohr = ^/{am), which is large compared 
to the Compton wavelength \q = 1/m and hence disparate scales appear allowing for an 
effective field theory treatment. However, the small scale 1/RBohr depends on the mass 
such that it does not approach a finite limit as m — oo, even for running as- 

For the heavy quarkonia we shall not assume any type of potential or binding model, 
but rather leave the matrix elements of the operators as non-perturbative entities, which 
may not be determined within the effective theory. The only thing that has to be kept in 
mind is that these matrix elements depend on the mass, as can be seen from the Coulombic 
example. The Coulombic system is, however, an extreme case, since all scales are set there 
by the mass of the constituents; we shall argue below that the mass dependence might, 
in reality, be much weaker. 

Based on the new equations of motion we get for tree-level matrix elements the rela- 
tions 

{tvD)hi+^ = ^^/iW = 0{l/m) , {tvD)hi-'^ = = 0{l/m) . (20) 

2m 2m 
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As we shall see, this means that the static equations of motion hold up to terms one order 
higher in the heavy-mass expansion. This is very similar to what happens in the v/c 
expansion of [Q, where the right-hand side of ( pOD is suppressed by an additional power 
in v/c. 

The magnetic- moment term of order 1/m appearing in the Lagrangian and all higher- 
order terms are treated as perturbations and lead to time-ordered products; these terms 
may be identified with the corrections to the states. The expansion of the fields Q^"^-* leads 
to corrections, which are in general local operators. However, as in any effective theory 



we are free to perform field redefinitions |Tl|], and hence the expansions (|g) and (y) are 



not uniquely defined; one may always move terms appearing in the expansion of the field 
into the Lagrangian; these contributions appear in the Lagrangian as terms that would 
vanish by a naive use of the equations of motion (|^). However, inserted into time-ordered 
products they yield local terms according to 

{^lj\T{U+^T{ivD)h''+^U+^T'h^+^}\'4j) = i6\x){4j\h'^+^rP+r'h^+^\ilj) + 0{l/m) , (21) 

where the higher-order terms appear, since the equation of motion is not quite the static 
one. However, this ambiguity of shifting contributions between the Lagrangian and the 
fields appears only in terms of order 1 / m? or higher and thus does not affect our arguments 
concerning the heavy-mass limit based on reparametrization invariance. 

Hence only the combination of (^ and (P) has physical significance and for the present 
application it is convenient to use expansions somewhat different from (^ and (|^). By a 
linear field redefinition we arrive at equivalent expansions, which are up to 1/m^: 

^ = -^static + (22) 



where 



-^static = ^^^^ ^^"^ - ^^"^ 



Here we have defined 





= 7^;+) 


+ Kt^ 




Gi 


= gS+) 


+ G[-^ 


g[^^ 






+ Kt^ 




G2 


= gw 


+ Gi~^ 


Gt' 



H)h^^^a^,{itD^), [i-tvD), (zZ}^)]}/.W . (24) 
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The corresponding expansion of the field Q\^^ reads 
QW(x) = (l + ^P^ilf - ^^{^vD)P.^I^ 



^ ( iiDf - iivDf - la.^iDHD'') + OillrrV^)] /i^ ix) . (25) 



In fact this is the form that has been obtained from QCD by a sequence of Foldy- 



Wouthuysen transformations . 

The advantage of this form is that now the Lagrangian no longer contains terms 
that would vanish by a naive use of the equation of motion; all these terms have been 
shifted into the expansion of the field Q\^^ ■ Again the corresponding expressions for the 
antiparticle fields are obtained by the replacement v —v. 

In this way one may obtain all terms which explicitly contain a heavy-quark field. 
However, there are more contributions appearing in order which are due to closed 

loops of heavy quarks in the full theory. These contributions may be expressed as local 
higher-dimensional operators involving only gluon fields. To order there are only 

two independent operators, leading at tree level to a contribution Cgiue to the Lagrangian 

( C.iue = TT{HD.,GnHD\GxJ} + -^^^ TriG^'^fG^^, J} (26) 

\2mJ ^ 307rm2 ^' ' ^^^^^ SGOnm^ ^ ^ ^' ^ ' 

where the gluon field strength is defined as 

[iD^,iD,]=igsG^, . (27) 

We may use the equations of motion of the gluon field to rewrite two of the 1/m? 
operators as 

h^+\{iD^), [{-ivD), (iD^)]]/i(+) = 47ra,/i(+)7^T"/iW 5] g7'^T'^g 

<? 

Ti{[iD,,G>^^][iD\Gy,]} = -2j2{qYT''q){q\Ty') . (28) 

qq' 

In this form they have a simple interpretation: The first one is the interaction of the 
heavy quark with the (virtual) light quarks in the quarkonium, and the second one is the 
interaction among these (virtual) light quarks, which is introduced by heavy-quark loops. 
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In order to calculate short-distance QCD corrections within this effective theory, one 
has to start from the Feynman rules as derived from Lq. The propagator H{k) obtained 
in this way is 

and hence contains all orders of 1/m. However, the reason why we had to include these 
higher-order terms was that this removes the divergent phase occurring in the limit of 
small relative velocity. This phase is a long-distance effect and may be absorbed into 
the states, which thus have to evolve according to the dynamics dictated by Cq given 
in ([T9|) . Once the phase is removed, we may expand the remaining expression for the 
short- distance contribution in powers of 1/m". Hence, as far as practical calculations are 
concerned, we may simply use the static propagator 

Hstat{k) = P+^^ (30) 
vk + le 

as in usual HQET, which is the leading term of (^). If we chose the velocities of the 
heavy quarks to be equal, then ill-defined imaginary parts such as 

"^'^ ^;5{vk)5{v'k) ^ / .^h5{vk)f asv-^v' (31) 



(27r)4A;2 ^ ^ ^ ^ J (27r)4 P 

will show up at the one-loop level, which are contributions to the divergent phase. How- 
ever, these are absorbed into the states and the real parts may be calculated simply in 
the static limit, even for equal velocities of the heavy quarks. 

In this way we may exploit the full machinery of the static limit, i.e. of HQET. However, 
although the short- distance contributions are calculable in the static limit, the matrix 
elements of operators composed of static fields will not be flavour-independent in the case 
of a quarkonium, since the states contain a non-trivial mass dependence. 

This concludes the set-up for the heavy-mass expansion for heavy quarkonia. In the 
following section we shall apply these ideas to inclusive annihilation decays of heavy 
quarkonia. 



3 Separation Between Long and Short Distances 

In this section we shall discuss the general set-up for a QCD-based calculation of the 
inclusive annihilation of a heavy-quarkonium state. The aim is to establish a separation 
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between long- and short-distance physics; the latter may be calculated in perturbation 
theory, while the long-distance part is non-perturbative and is parametrized in terms of 
hadronic matrix elements. 

The starting point is the inclusive transition rate F of a quarkonium state I?/'), which 
is given by the optical theorem in terms of the forward matrix element of the transition 
operator T 

r = 2 Im (V'|T|V') , (32) 

which is itself related to the discontinuities across the cuts of the two-point function of 
two fields interpolating the state I?/') 

G{p) = J ^e-^^^(0|r{^(x)^(0)}|0) . (33) 

For the case at hand, (|33D has a cut along the real axis of the complex p'^ plane starting 
at mf^gf^^, where the mass mught corresponds to the lightest hadronic state the heavy 
quarkonium may decay into. 

On the other hand, the transition rate F is related to the discontinuity across the 
cut at the mass of the heavy quarkonium state, which is much larger than mught- If we 
consider only decays into light hadrons we are thus far away from the resonance region, 
for which is of order mf^gf^^. Using local duality we may calculate the annihilation of the 
two heavy quarks into light hadrons in perturbation theory, since the scale for this process 
is set by the heavy-quark mass. In other words, the short- distance piece of the transition 
rate F has an expansion in inverse powers of the heavy-quark mass, up to logarithms 
induced by renormalization, and the coefficients of this expansion may be calculated in 
perturbation theory. 

The annihilation part of the transition operator T in terms of the heavy-quark fields 
is in general given by 

T=[d'Xd'pd\ Y.Y.K\f\x,p,i)Q{X + p)VfiQ{X-p)Q{-i)VuCQ{i) , (34) 

C=l,8 j,k 

where C is a matrix in colour space; C = 1 (C = 8) corresponds to the colour combination 
C0C=101 (C0C = T"{g) T") and the sum over j and k runs over the sixteen Dirac 
matrices. The kernel Kjk depends on the c.m.s.-coordinate X of the heavy-quark pair as 
well as on the two relative coordinates ^ and p. 

To set up an expansion in inverse powers of the heavy-quark mass we perform a phase 
redefinition of the heavy-quark fields 

Q{x) = ex.Y){—im{vx))Q[^\x), Q{x) = ex.]i{—im{vx))Q[^\x) (35) 
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where the superscript (+) ((— )) refers to the annihilation part of the fields for quarks 
(antiquarks). Here v is a velocity vector, which we chose to be the same for both the 
quark and the antiquark. Hence we shall work in the limit we have discussed in the last 
section and finally identify the velocity with the velocity of the heavy quarkonium. 

This phase redefinition removes the dominant piece of the space-time dependence of 
the heavy-quark field operator; the remaining dependence is only due to the residual 
momentum of the heavy quark in the heavy hadron. 

Inserting these redefined fields the transition operator takes the form 

•' C=l,8 j,k 

qW(x + p)TfiQ[r\x - p)g(-)(-orfcCgW(0 . (36) 

Introducing the Fourier-transform of the kernel as 

(X, M) ^ / |^^|^e-^-e-e.-K<?(P, ,) , (37) 

we obtain 

T . j,'X,W% S L/ ||^(S,|^^5?'(A-..)e-«<^--'e--e- 
Q^+>(X + p)TjCQl\X - p)(3<->(-0rtCQ(+>K) . (38) 

The dominant momentum dependence is 2mv and we can remove this large piece by 
redefining the c.m.s. momentum; hence we introduce the residual c.m.s. momentum 

P = P - 2mv (39) 

and obtain 



C=l,8 



gW(x + p)TfiQ[-\x - p)gl-)(-Of fcCgW(o , (4o) 

where the kernel k!"^ now depends only on small, i.e. residual, momenta and on the heavy 
mass m. 
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Furthermore, the kernel has an expansion in inverse powers of the heavy-quark mass m, 
up to logarithms of the form ln(m//i), where /i is a factorization scale0. The transition rate 
r is independent of the factorization scale, and thus the scale dependence of the expansion 
of the kernel has to be compensated by the corresponding dependence of the matrix 
elements, which appear in the expansion. As usual, this scale dependence is governed by 
the renormalization group, and hence one may recover the logarithmic dependence of the 
transition rate on the heavy-quark mass m. 

The expansion of the kernel thus takes the general form (A = fi/m) 

lCf,\m;P,n,v) = <^(A) + ^ (P,sf ^(A) + vr,T,f ^(A) + r/,f/f ^(A) 



2m 

+ (^)' {P.P.Slf '''''' (X) + P.TT.sf (^)-(A) + P,,.5f (^)-(A) 

+ 0{l/m^) . (41) 

The powers of the momenta will simply yield matrix elements of operators involving 
derivatives. For later use it is convenient to introduce those that correspond to the 
residual-c.m.s.-momentum (RCM) 



^^, (gl+^rg(-)) = (gWr(^D,)g(-)) + [Qi-^\^ ^,)rg(-)j (42) 

and to the residual-relative-momentum (RRM) 

'Q^'\^ ^,)TQi-A = (gl+)r(zZ}^)g(-)) - fgl+H^ D,)TQi-A (43) 



where D = d — igA is the QCD covariant derivative for fields in the fundamental rep- 
resentation. For the other bilinear gi~'* ' ' ' Qi^^ we define the derivatives in the same 
way. 



^We restrict our attention here to the purely perturbative aspects of the operator expansion and ignore 
possible problems induced by renormalon poles. These have been discussed recently in the context of the 
heavy mass expansion in [fsl. 
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Inserting (^) into the expression for the T operator and using the definitions (|42| ) and 
^) one finds for the transition rate F up to and including order Xjvr?: 

m^)= E E|^5?W^)(^l[01+^r,cQ(-)][Q(-)f,cQW]|^), 

+ ^ [5f '^(A)(V^|[z9,(QWr,CQ(-))][gi-)f,C'QW]|V'), 

+ T,f-(A)(^|[(QWr,(^ D^,)cgi-)][Q(-)f,CQW]|^), 
+ f^f "(A)(^|[(Q^)r,CQ(-))][Ql-)r,c(^ d^,)qw 

/ 1 \ 

+ 



+ (^)'^^(A)(v^|[^9,(QWr,(^ :d .)cg(-))][g(-)r,cgw 
+ 5f (^)^^(A)(v^|[za,(gWr,cg(-))][g(-)r,c(z D^.)gW]|^) 
+ T,f (^)'^'^(A)(v^|[gWr,(z :D^,)cg(-))][z9.(g^)r,cgw 
+ T,f (^)^'^(A)(^|[gWr,(z :d^,)(^ :D^.)cg(-)][g(-)r,cgw 
+ r,f (^)^'^(A)(v^|[gWr,(z :D^,)cg(-)][g(-)f,c(^ D^.)gw 
+ f7,f (^)''^(A)(^|[z9.(g(+)r,cg(-))][gl-)r,c(z D^,)gl+) 
+ f7,f (^)^^(A)(^|[gWr,(^ :^.)cg(r)][gl-)f,c(^ D^,)gw 

+ t/,f ^'^^^(A)(^|[gWr,cg(-)][gl-)(^ ^.)(^ 



+ C(l/m3) (44) 



where the subscript /i at the matrix elements indicates their renormalization point and 
all field operators have to be taken at x = 0. In what follows we shall not display 
the dependence on x any more, if x = 0. Furthermore, we have replaced the ordinary 
derivatives appearing in the Taylor expansion of the matrix elements by covariant ones in 
order to ensure gauge invariance; this could be implemented from the very beginning by 
defining the T operator including appropriate Wilson-line operators. 

The lengthy expression ( ^4|) is the most general short-distance expansion for the matrix 
elements of the T operator for the case of heavy quarkonia decay up to and this 

is the first ingredient for a 1/m expansion of the transition rate F. In a second step, one 
has to expand the matrix elements appearing in p^ ) in powers of 1/m. 
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4 1/m Expansion for Inclusive Quarkonia Annihila- 
tion 



In this section we write down the complete heavy-mass expansion for the matrix elements 



appearing in (0). The result we are aiming at is an expansion of the transition rate T of 
the form 

r = ro+(i)r. + (i/r. + -- (45) 

where we shall explicitly construct the terms up to and including In this expansion, 

the ratio of two successive Tj is set by a small scale 



A (46) 



where A is small compared to the heavy-quark mass. However, there is some difference 
between the heavy quarkonia and the heavy-light systems. In the latter the parameter 
corresponding to A is A = Mhadron — which becomes independent of the heavy quark 
mass in the heavy-mass limit. 

As we have discussed above, in heavy quarkonia systems this is not true, since in order 
to perform the limit v w we had to include subleading terms of the 1/m expansion into 
Cq, which determines the dynamics of the states. In this way a non-trivial dependence 
on m is introduced into the scale A induced by the binding of the two heavy quarks. 
Nevertheless, this is a small scale and the expansion will be useful. 

Since this means that we are not able to perform a static limit for a heavy quarkonium, 
we shall no longer have heavy-flavour symmetry. However, if we define the scale A in a 
similar way as A for heavy light systems 



zvd{0\hi-W+^\ij{v)) 

(o|/il-^r/il+VM) 



we may get some idea of the flavour dependence by looking at the level spacings of the 
(n^Si) quarkonia of the c and the b quarks. Here it is a well-known fact that the level 
spacings are almost the same in the (cc) and [bb) system, which we take as a hint that 
A is only weakly dependent on the heavy flavour. We shall return to this point when we 
discuss heavy-quark symmetries in section 7. 

The way we shall set up the expansion is to use the expansions (EBf) and 



The fields Q^,^-* are expanded as in (UB) while the chromomagnetic term Gi and all terms 
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of order 1/m? of the Lagrangian are treated as perturbations, which lead to time-ordered 
products. Up to and including l/rn? only the double insertion of Gi and the single 
insertion of the order 1/m? terms of the Lagrangian will play a role. 

The leading-order term Fq in the heavy-mass expansion is given by forward-matrix- 
elements of dimension-six operators involving the static fields h'^^\ In general, there are 
four dimension-six operators, two with colour structure 1 {S> 1 



= {h^+h^.h^'^) {h^-h^h^^^) (4J 



and the corresponding two operators A\ , A2 with colour structure T'^^T'^. In fact, 
due to the projections P± = (1 ± ^)/2 implicit in the static operators these are the only 
spinor structures that can emerge (see appendix for details). 

Thus the leading term Tq is always given as a linear combination of the four dimension- 
six operators A] (the factor 1/2 arises from (p2D) 



oro = E E c{AP){mP\i^) , (49) 

^ i=l C=l,8 

where the coefficients C{Af^^) depend on the specific process, i.e. the kernel Kij appearing 



in (|4^) . We shall give a few examples of this "matching procedure" below. 

In higher orders of the heavy-mass expansion matrix elements of higher-dimensional 
operators appear from the expansion of the , which have increasing powers of covariant 
derivatives acting on the static fields and it is convenient to use the RCM and RRM 
derivatives introduced in (^2)) and (^). 

All dimension-seven operators either vanish due to symmetries or they are proportional 
to (ivD). The latter may be rewritten using the equations of motion; in the static case 
their tree-level matrix element vanishes, while one obtains a dimension-eight contribution, 
if we include the higher-order terms of the equation of motion (^. Whatever choice one 
prefers, there are no tree-level matrix elements of local dimension-seven operators. The 
only contribution which is allowed at order 1/m is a single insertion of the first-order 
chromomagnetic-moment operator 



|t{[gS- 



'^^'''x)+G[-\x)[aP}\^) (50) 



for which we shall show below that it vanishes due to spin symmetry. Thus for the 
inclusive heavy quarkonia decays there is no contribution of order A/m. 
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The second-order term in the 1/m expansion has contributions from both the La- 
grangian and the local operators. The local dimension-eight operators may be classified 
according to the derivatives defined above. The first set of operators consists of the contri- 
butions with two RCM derivatives (RCM x RCM). There are in general six independent 
contributions of this type, three of which are colour 101 and the other four are the 
corresponding T" operators. The colour 101 operators are 



(1) 



B. 



(1) 

2 

(1) 



id^, [U^'^Yh^''^)] [id'' {U-^-i,h^+^ 



(51) 



Note that one may flip the RCM derivatives from one heavy-quark bilinear to the other, 
thereby picking up a total derivative, which will not contribute to the forward matrix ele- 
ments we shall consider. This is the reason why no terms appear in which both derivatives 
act on the same heavy-quark bilinear; these may be rewritten into the B operators. 
In the same way one obtains six colour 101 RCM x RRM operators 



a 



(1 



6 — 



^(+^5(^ Dt,)h^-^j [id" {h^-^^5h^+^)] +h.c. 

V)^^"^) [id" {U-h^^h^^'^)] +h.c 

+h.c 
+ h.c 

7/3(^ ^.)h^-^) [td,{h^-h5h^+^) 



IS' 



)d, {h^-hph^^^) 



-Fh.c. 
+ h.c.| 



(52) 



and the corresponding T" T" operators. 

Finally, there are 22 RRM x RRM operators, half of which are 101 and the other half 
are the corresponding colour T"0T'". These operators fall into two categories, namely 
(for the colour singlets) 



D^' = ^,)h^-') {h'-h.{^ ^')h^^') 
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D 



D 



D 



(1) 



(1) 



(1) 



and 



(1) 



(1) 



(1) 



(1) 



h^^h^ii 7)u)h^-^^ Z^V^^^ 
h^^huit iT^)/^^"^) (h^'h^i^ l^'^)h^^^) 
le'^f'^^v^ (h^-^hp{i Zr«)/i(-)) (h^-h^i^ lD,)h^^^ 



h.c, 



(53) 



h^^h^ii D fh'^'^j [U-^^5h^+'>) + h.c. 



(54) 



Thus there are in general 40 dimension-eight operators. Note that we have aheady 
dropped the operators which vanish because of the equations of motion for the heavy 
quark. 

The general expression for r2 is thus given by a linear combination of the 40 local 
dimension-eight operators and the non-local contributions from the time-ordered products 
with the Lagrangian (recall the factor l/(4m^) in the definitions (|^), (P3|), and (^) and 
the factor 1/2 in (|3|)) 



6 



i=l C=1,S i=l C=l,8 

+ E E c{DP){ij\DP\ij) + Y: E c{EP){m 

i=l C=l,8 i=l C=1,S 

+ E E CiAP)i-t) ld'x{^\T{L,ix)Af\0)}\^p) 



=1 C=l 
\2 2 



+ 



1=1 C=l,i 



E E C{aP) / d'xd''y{ij\T{G,ix)GMA. 
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(a) (b) 

Figure 1: Fermion loop diagrams. Diagram (a) contributes to the electromagnetic kernel 
for the quarkonia decays into e~^e~, diagram (b) contributes to the strong interaction 
kernel for the quarkonia decays into light hadrons. The dashed line means that the 
intermediate states are on-shell. 



+ E E C{aP){~z) / d'x{^\T{C,Ux)AP{0)} \^) (55) 

i=l C=l,8 

where the coefficients in the time-ordered product terms are again given by the lowest- 
order coefficients C{aP). This expression is still quite general and simplifies once the 
heavy-quark symmetry and the fact that we have to take forward matrix elements is taken 
into account. 



5 The Operator Coefficients to Leading Non- Trivial 
Order 

In this section we shall calculate two simple kernels in order to show how our method is 
applied, namely the fermion loop diagrams depicted in fig. |^. 

These diagrams yield contributions of order a (fig. |l|a) and as{m) (fig. |I]b). Calculating 
the electromagnetic contribution (fig. ^) one obtains for the kernel 

4:71 

/C^^(P,7r,r7) = -^a^c^g^u , (56) 

where cq is the charge of the heavy quark in units of the electron charge. Furthermore, 
P = p + p is the sum of the momenta of the quark p and the antiquark p. Hence only the 



Dirac matrix combination 7^ ® 7^ contributes in the sum over the Dirac matrices in (@). 
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The kernel is expanded using 




and the result is matched to the operators by identifying the momenta with the deriva- 
tives appearing in the operators (see appendix for more details). The only non- vanishing 
coefficients are 





m) 




Sm^ ^ 






m) 




^ 2 2 






m) 


- ^ a^c^ 




(57) 



From this one may obtain the coefficients for the quark-loop diagram of fig. |I]b by the 
replacements Cga^ — nja^(m)/2 in the coefficients and 1 (8) 1 ^ in the operators, 

where is the number of light flavours which are allowed in the quark loop. In this way 
we obtain 



c««(4^\ 


m) 


= C''i{Bf,m) = 






C«^(Ef\ 


m) 


= C''\Ef\ni) = 






C«^(Ef\ 


m) 


= 12m2"'^"^^^^ 




(5^ 



The matching calculation yields these coefficients at the scale n = m and hence it is a,^ 
taken at this scale that enters the expression. We have indicated the /i dependence of 
the coefficients by an additional argument for these functions. A change of this scale is 
governed by the renormalization group of the effective theory, which is discussed in the 
next section. 



6 QCD Evolution of the Coefficients C 

In general QCD corrections render the operators and the coefficients scale-dependent in 
such a way that the transition rate F is scale-independent. Schematically this may be 
written as 

T = Y.C^i^^)ml (59) 

i 

20 



where ()|^ means that the matrix element is normahzed at scale /i. The scale dependence 
of the matrix elements and the coefficients Ci is governed by the anomalous-dimension 
matrix 7 

fx-^mi = -^,,{o,)i, (60) 

which is obtained in the standard fashion from the ultraviolet divergences of the matrix 
elements. In order to have the transition rate scale-independent, the coefficients Cj(yu) 
have to obey the renormalization group equation 



(61) 



where the initial condition at the scale fi = m is given from the matching calculation, i.e. 
the calculation of the hard kernels as performed in the last section. 

(C) 

The four dimension-six operators A\ have a vanishing anomalous- dimension matrix 
and hence there is no mixing between these four operators. Keeping the velocity of the 
quark operator different from the one for the antiquark we find that the anomalous- 
dimension matrix - after the redefinition of the states - vanishes as 



(6) / \ 

7ij oc ujr{uj) 

where lu is the product of the two velocities and 

1 



r{x) 



ln(x + — 1) . 



(62) 



(63) 



Hence the coefficient functions C{Af^'') are scale-independent. 

In a similar way we find that there is no mixing between the local dimension-eight 
operators. In order to access the renormalization of the non-local time-ordered product 
terms, we have to ffist study the renormalization of the terms in the Lagrangian. The 



renormalization of the ffist-order Lagrangian is [T^ 



C0C2 G] 



C2 = ry-9/(33-2n^) 



where 



as{m) 



(64) 



(65) 
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(a) (b) (c) 




Figure 2: Feynman diagrams leading to mixing between the non-local and the local 
operators. The square boxes correspond to insertions of the first-order kinetic energy 
operator Ki and the shaded blob corresponds to the dimension-six operators A] . The 
diagrams (a) to (e) have to be calculated for all permutations of the external lines. 

and Co is the wave function renormalization constant of the static field, which reads in 
Feynman gauge 

The subscripts n and m label at which scale the matrix elements, in which the operators 
are inserted, have to be normalized. We note that Ci = 1 ensures that both terms 
of the Lagrangian (|T^) renormalize in the same way; in other words, the wave-function 
renormalization is the same in both the static and the effective theory based on (^9|). 

Of the second-order Lagrangian we need only the contributions -^'2^'' and Cgi^e, since 
there will be no contribution from a single insertion of due to spin symmetry. Fur- 
thermore, we shall not include the renormalization of the second-order Lagrangian, since 
its contribution is hard to estimate, no matter at which scale we consider the matrix 
elements. 

One important effect of the renormalization in the present application is the mixing 
of the time-ordered products of the Lagrangian into the local operators, since non-local 



22 



operators may require a local counterterm. The only contribution that leads to an ultra- 
violet divergence is a double insertion of the first-order kinetic energy operator Ki. We 
define the non-local operators 



T 



d^xd^yT [Ki{x)Ki{y)Af\0)} 



(67) 



and the one-loop diagram types leading to ultraviolet-divergent contributions are depicted 
in fig. 0. Diagrams (a) to (c) are one-particle irreducible pieces, while [d) and (e) may be 
interpreted as the mixing of the time-ordered product T{Ki{x)Ki{0)} into local terms of 
the Lagrangian. 

We define two sets of seven operators, the first corresponds to the spin singlet, the 
second to the spin triplet coupling 



{Ti'\T}'\D['\E['\B['\D^^\Ei''> 



(1) T^») riW r(8) r)(s) 



(6^ 



These two sets have the same anomalous-dimension matrix and we find 



67r2 












-2 





4 











8/9 





-6 


5/3 


-7/6 












































































































(69) 



Using the anomalous- dimension matrix, one may calculate the scale dependence of the 
coefficient functions 



(cfi) = (c(Tr),c(Tr),c(Dro,c(i?r),c(sr),c(z^r),c(£;r)) 

(Cfl) = (C(T«),C(Tf ),C(Z}(^^C(E«),C(i?f ),C(Df),C(Ef )) 



.(l)^ 



1(1) 



?(8)^ 



1(8) ^ 



^(8). 



(70) 



by solving the renormalization group equation (^) with the initial conditions given at 
/io = m. The solutions may be expressed in terms of the value of the coefficients at the 
scale m and one obtains 

C{Ti'\^i) = C{T^\m)=C{A'l\m) 
CiTf\f^) = CiTi'\m)=CiAf\m) 
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C{D?,^^) = C{D^^\m) + f ^^^^^(Tf , m) Inr^ 

66 — ZUf 

C{Bf\fi) = C{Bf\m) - —^^C{Ti^\m)\nr] 



33 


-2nf 




24 


33 


-2nf 
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C{Df\^,) = C{Df\m) + — -^C(T«,m)lnr^ 

33 — in/ 

H C(t|^'', m) ln?7 

3 33 - 2n/ ^ ^ ' ^ ' 

C(Ef\/i) = C{Ef\m) - y^^^C(T(^),m)lnr/ (71) 

and the same solution for the second set of coefficients. 

The non-trivial evolution of all the coefficients of the local operators is driven by the 
coefficient of the time-ordered product corresponding to a double insertion of Ki\ in turn, 
this coefficient is simply the one of the dimension-six operator with the corresponding 
spin structure. 

The scale /xq is the matching scale, which we chose to be /iq = m. On the other 
hand, we want to study the matrix elements of the operators at some lower scale, say 
the scale of the binding energy of the quarkonium, where we expect to have a reasonable 
approximation using a wave-function model. The evolution equation ( |6T1) allows us to 
change the renormalization scale, thereby inducing mixing of the non-local operators into 
the local ones. 



7 Forward Matrix Elements and Heavy Quark Sym- 
metry 

Once we have scaled down to the small scale of the order of the binding energy of the 



quarkonium, we have to evaluate the forward matrix elements appearing in (55). This 
requires in general non-perturbative input, which has to be supplied by other methods 
such as lattice gauge theory or by model estimates. However, we may use heavy-quark 
symmetry to count the number of independent non-perturbative parameters. 

The Lagrangian given in (|19|) still has the heavy-quark spin symmetry. As a 
consequence, the heavy quarkonia systems fall in general into spin symmetry quartets: 
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For a given orbital angular momentum £ and radial excitation quantum number n, the 
four states (in the spectroscopic notation '^^^^£j) 

[nH, n^£e-i nH, nH,+^] (72) 

form such a spin symmetry quartet. An exception are the S waves (£ = 0), for which the 
three polarization directions of the n^Si and the u^Sq form the spin symmetry quartet. 

In order to exploit the consequences of the spin symmetry for the inclusive hadronic 
decays we shall use the trace formalism. We denote with \Yi) the spin symmetry quartet 
consisting of the spin singlet and the spin triplet for a given orbital angular momentum 
£. The coupling of the heavy-quark spins may be represented by the matrices 

, , J P+75 for the spin singlet , . 

^Y{V) j for the spin triplet ' ^ ' 

where P+ = (l + |^)/2 is the projection of the "large components". Note that the matrices 
Hy are independent of i. 

For the dimension-six operators this implies 

{Ye\ (h^+^CTh^-^) [U-^CT'h^+^) \Ye) = af^ Ti{HyT] Ti {V Hy} , (74) 

which means that there are only two independent parameters for a given £ describing 
the dimension-six matrix elements, namely a^p and In the study of exclusive non- 
leptonic decays of heavy-light mesons, matrix elements of four-quark operators are usually 
estimated in vacuum insertion. Applying this approximation to (|7^ suggests that a^^ is 
the dominant coefficient. 

The generic dimension-eight colour 101 operators may be written in terms of these 
representations as 

{Yi\[id^ (/iWr/i(-))][i9" (/i(-)r'/iW)]|F^) = fef Tr{;^yr} Ti{V'Hy} 
{Yi\ {h^+\i ir'')r/i(-)^ [id"" (/i(-)r'/iW)]|F^) + h.c. = cf Ti{HyT] Ti{V'Hy} 



{Y,\ {h^+\i ir'')r/i(-)j (h^-^r{i lD'^)h^+^J]\Y,) = Tr {ifyF} TrjF'i/y} 
{Y,\ (U+\t 7)''){t 7)'')Th^-A {h^-^T'h^+^)]\Yi) + h.c] = ef Tr {HyT} Tt{T'Hy} . 



(75) 

The tensors b'^'^ , c^'^, d'^'^ and e^'^ have to be constructed from the velocity vector v and the 
metric tensor. Furthermore, if one contracts one of the indices with v, one may use the 
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equation of motion of the fields /i^^^; in the static hmit this simply vanishes, while with 
an equation of motion as ( p!8|) one obtains a term one order higher in the 1/m expansion, 
which may be dropped, since we work only to order 1/m?. Consequently, the tensors are 
given by 

= dfXg^''' - v^'v") , = eeig^"" - v^v"") . (76) 

Spin symmetry thus places a very strong restriction on the local matrix elements, since 
for a given orbital angular momentum i there are only eight non-perturbative parame- 
ters 6^ , ■ ■ ■ e^l^ and the corresponding octet partners bf^ , • • • ef^ , which describe all the 
forward matrix elements of the local dimension-eight operators given above. 

From an estimate based on vacuum insertion one would guess that fog^'*, and Cq^^ 

(C) 

dominate all other coefficients; in particular all the expected to be small, if vacuum 

insertion makes any sense. 

In a very similar way one may analyse the spin symmetry structure of the non-local 
terms. The kinetic energy terms are spin symmetric; furthermore there is no contribu- 
tion from insertions of Ki, since this is already contained in the dynamics of the states; 
however, we note that this does not mean that Ki does not show up at all: A double 
insertion of Ki mixes under renormalization into local dimension-eight operators, but this 
is a short-distance effect which was calculated in the last section. 

For an insertion of the second-order kinetic energy term K2 we obtain from the trace 
formalism 



—I 



J d^x {Ye\T{K2{x) (/i(+)r/i(-)) (/i(-)r'/i(+)) \Ye) = k'i^ Tr [HyT] Tr {VHy] (77) 



and a corresponding expression for the colour T*^ T" contribution. 

Furthermore, spin symmetry implies that a single insertion of a chromomagnetic mo- 
ment operator vanishes. We have {j = 1, 2): 



d^x {Y^\T{Gj{x) (/i^r/i^-)) (/i(-)r'/iW) \Ye) = 

Tr {HYa^,P+T} Ti{T'Hy}+ Tr {TlYTP^a^,} Tt{T'Hy} 
^Tt{HyT} TT{T'P+a^,HY}+ Tr {HyT} Tr {a^.P.F'ify}] . (7J 



Thus the tensor g^'^ has to be antisymmetric, but it may only be built from the velocity 
vector and the metric tensor; from these only symmetric combinations are possible, and 
hence these contributions have to vanish. 
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Finally, there is also a double insertion of the chromomagnetic moment operator Gi. 
Analysing the spin structure of this term gives 

^ J d'xd%{Y,\T{G,{x)G^{y) (U+^Th^'^) (/i^^^r'/iW) \Ye) = S^'"'^ 

' Tr {Hyo^^P+g^pP+t] Tr {T'Hy] + Tr [HYP+a^^VP^a^p] Tr {T'Hy] + ■ ■ \n) 

where the ellipses stand for all other possible insertions of the two sigma-matrices. The 
tensor S now has to be antisymmetric and it may only be built from the metric tensor 
Qfj^u^ since we have P±a^i,v'^P± = as well as P^a^j^yV^P^ = 0. Hence we have only a single 
parameter for this contribution 

J d'xd'y{Ye\T{G^ix)G^iy) {h^+W-^) {h^-^T'h^^^) = G« 
"Tr {HYa^yP+a^'-P+r} Tr {r'ify} + Tr {:^yP+a^,rP_a^^} Tr {r'ffy} + ■ ■ •](80) 
and a corresponding relation for the colour combination T'^^T"'. 

(C) (C) 

The constraints from spin symmetry thus allow us to drop the operators C5 , Cg , 

(C) (C) (C) 

in a calculation for the total rate. Of course, the contribution of Cgiue, 
(PB|), is spin-symmetric, since it involves only the light degrees of freedom. Furthermore, 
the only contribution of chromomagnetic-moment operators is the double insertion of Gi, 
which in turn is the only spin-symmetry-violating contribution that appears. 

At the end of this section we want to discuss the flavour symmetry. Let us recall 
that for heavy-light systems both the spin and the flavour symmetries hold. As already 
mentioned several times, the flavour symmetry is broken in the present case of heavy 
quarkonia, because the states acquired a mass dependence. The interesting question 
arises about the amount of flavour-symmetry violation. As a simple example, we study 
the simpler problem of flavour-symmetry breaking for mass splittings. In the heavy-mass 
limit the four states in ( |72D should be degenerate; splitting between the members of 
the spin-symmetry quartet is induced by spin-orbit and spin-spin interactions. We may 
analyse this by looking at the mass of a heavy quarkonium in terms of the Lagrangian 

+H) (^)7 d'x {^1^\T{[g[-'\x) + g[~\x)][G[^^ + G[-^]m 
+0{l/m^) . (81) 
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Note that no insertions of k['^^ appear, since this term is already included in Cq. 

The first contribution A{n,i) is the "binding energy", as one would obtain from the 
solution of the Schrodinger-type of equations corresponding to Cq. This parameter can 
only depend on the quantum numbers n and i, since spin symmetry is still unbroken. 

Breaking of spin symmetry occurs first at order 1/m; it is given by the expectation 
value of the first-order chromo magnetic moment operator. In a non-relativistic language, 
this would correspond to a coupling of the form (s^^^ + s^~^) ■ B, where B is the chromo- 
magnetic field, which is created by the orbital motion. Hence we expect B cc L, where 
L is the orbital angular momentum. In other words, the first-order term is a spin-orbit 
coupling term. 

The spin symmetric local terms of order are corrections to A{n,i), while the 

second-order chromomagnetic operator yields a correction to the spin-orbit term. The 
time-ordered product of the two first-order chromomagnetic moment operators will give 
(aside from a correction term to A{n,£)) a spin-spin coupling term, which in a non- 
relativistic form is s^^^ ■ s*-^-* for £ = 0. Hence we arrive at a mass formula of the form (for 



0): 



M(n2^+^£j) = 2m + A(n, i) + ^4?^ I [J{J + 1) - ^(^ + 1) - S{S + 1)] 

2m 2 



r(n,. 



(82) 



where uj{n, €) and r(n, t) are given in terms of matrix elements involving the Lagrangian. 

Spin-dependent relativistic corrections can be systematically studied in the Wegner- 
Wilson loop approach. Their general structure up to order (under the assumption 

of a fixed background field A^{x), cf. ref. |T^) has been obtained by Eichten and Feinberg 
(EF) 

... ^ , , V,{ry + 2V,{r)' + 2V,{ry ^ ^ V,{r) ^ ^ 
V{r) = V,{r) + L.S + ^s,.s, 

V^{r) (f- Sir -32 1 ^ ^ \ 

+ ~ ~2 o ^1 ■ ^2 

-\- spin-independent corrections , (83) 

where Vo(r) is the spin-independent potential and Viir) are related to expectation values 
of the colour-electric and magnetic fields, e.g. (r = ri — r2 and the limit T — > oo is 
understood): 

(t?-^) ^3(r) + ^^4(r) = |J f_l'^ f_''^^ (i?.(ri,t)5,(rl,t'))dtdt'/(l) ■ (84) 
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While (^) is exact for the spin-dependent relativistic corrections (for equal masses), the 
question of the spin-independent relativistic corrections appears not to be yet completely 
understood at the present time (for a review see e.g. [0). 

The potentials Vi{r) are usually assumed to be flavour-independent. This assumption 
is not only in agreement with observed mass splittings of charmonium and bottomonium, 
but also confirmed by lattice-QCD studies (for a recent study see e.g. ref. fl^l)- The flavour 
dependence of mass splittings in the EF formulation (^) enters through the explicit 
factors, but also through the mass dependence of the expectation values 

(r"? Vl''\r))ni = r dr RUr) r« ^/^^(r) , (85) 

JO 

where Rni{r) are the radial wave functions (normnalized to /o°° dr Rnii^) = !)• Our 
result (|8l| , |82D has a structure similar to that of (pSf ). For example, we may identify 

r{n,l) = {lv,{r))^i . (86) 

On the other hand, it is straightforward to include one-loop corrections to our result 
(pl|), while the EF representations seem to be non-renormalizable objects, and hence any 
attempt to calculate them in one- loop order should end up with a divergent result [p!5|] . 

Although the functions a;(n, £) and r(n, i) are matrix elements of flavour-independent 
operators, they in fact depend on the flavour through the states. We may obtain some idea 
on the flavour dependence of u!{n, i) and r(n, i) by comparing the spectra of bottomonium 
and charmonium. For example, the spin-orbit term (proportional to L ■ S" in (|83D ) and the 
tensor interaction (last term in (|83D ) do not contribute to either of the mass differences 
m{2^Si) — m{l^Si) or rh{l^Pj) — m(l^S'i), where fh{l^Pj) denotes the centre of gravity 
of the three ^Pj {J = 0, 1,2) states. Experimentally these mass differences are 563 MeV 
(589 MeV) and 440 MeV (428 MeV), respectively, for bottomonium (charmonium). Hence, 
certainly A{n, I) and T{n, I) are basically flavour-independent. 

Thus the parameters entering (^) seem to be only weakly dependent on the heavy 
quark flavour for the observable states of heavy quarkonia, although heavy-flavour sym- 
metry is not present any more once we include subleading terms into Cq. However, it 
is generally believed that in the heavy-mass limit a quarkonium should behave as an al- 
most Coulombic system. This would imply that A(n) ~ m, uj{n) ~ and r(n) ~ m^, 
which is not compatible with the data from charmonium and bottomonium. (For (|83|) , 
the vector- Coulomb case means spin-independent terms are proportional to a^m, while 
all spin-dependent terms are oc ajm.) This could mean that either these systems are not 
close enough to the heavy-mass limit to become Coulombic, or that the heavy-mass limit 
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is not the Coulombic one. Let us recall that the heavy-light systems with a b quark and 
to a lesser extent the ones with a c quark already behave as one would expect in the 
heavy- mass limit. This clearly indicates that the flavour dependence of the matrix ele- 
ments parametrizing the long-distance effects is smaller than suggested by the Coulombic 
limit. 



8 Some Phenomenology 

In this section we shall consider a few applications of the formalism outlined above. We 
shall not give a full discussion of the phenomenological applications of the heavy-mass 
expansion for heavy quarkonia, but rather study a few simple examples in some detail, 
based on the matching calculation performed in section 5. 

Compared to the conventional approach, and also compared to the v/c expansion 
advocated by BBL , we have a much larger list of parameters, although spin symmetry 
reduces this number to some extent. In a given order of the 1/m expansion they are all 
of the same dimension and there is no a priori reason why some of them should be less 
important than some others; at least the 1/m expansion does not give any hint. 

The most simple example (although maybe academic) is the semi-inclusive decays of 
a heavy quarkonium into e~^e~ and light hadrons. In leading order of the expansion we 
have that rjQ — >■ e~^e~ light hadrons vanishes, while we have for the ipQ 

^Toi^PQ ^ e+e- light hadrons) = C^^(4'\ /i)(V^l4'^l^) = ^«'4«o'^ (§7) 

with 

C-iA^\^^) = -^ycl. (88) 

There are no corrections of order 1/m, since there are no local dimension-seven oper- 
ators having non-vanishing matrix elements, and also the only non-local term, which is 
an insertion of a first-order chromomagnetic moment operator G^'' , vanishes due to spin 
symmetry. 

The first sub leading contribution appears at order 1/m^. Using the matching cal- 
culation of section 5 and taking into account the running considered in section 6, we 
obtain 

It, = C-(4^\^)(^|4^)|^), + C-(Df ,/i)(^|Df I^), 

+C-(E«,/i)(^|E«|^), + C-(E«,/x)(V^|E«|^), 
-Fr^(E«,/i)(^|EfV). + C^^(4'\/i)(A:?^ + i^^')) (89) 
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where /i is the renormahzation point /i < m, and F^^'^ is the contribution of the single 
insertion of the purely gluonic piece Cgiue of the Lagrangian. 
The coefficients are given by 



C^"(Df ) 
C-(4'\/i) 



71" 2 2 
TT 



^'4 

2 2 16 



ln?7 

2 2 

^22! 

tO; Cn 

3m2 ^ [2 33 - 



In?? 



(90) 



Note that we have an additional factor l/(2m)2 in front of the second-order contribution 
according to our definition p5|). 

This expression looks relatively complicated, but it simplifies somewhat due to spin 
symmetry. The matrix elements of the local operators are given in terms of the three pa- 



rameters 



el^'' and 



Q , while the non-local terms introduce another three parameters: 
2 



k'2' ■, G'^^\ and At order l/m^ there are thus six parameters describing the decay 

rate, which may then be expressed as 



1 



IT 



.,2„2 



126(^) + 8e« 



12 



+124'^ 



16 



33 - 2nf 



In 7] 



33 - 2n/ 



(91) 



A further simplification can only be achieved with additional theoretical prejudices. 
A popular assumption (although ad hoc) for non-leptonic inclusive decays of heavy-light 
mesons is vacuum insertion. Note, however, that this is not a scale-invariant statement. 
Usually vacuum insertion is applied at some small scale, where the matrix elements are 
estimated by e.g. a wave-function model. There the local contributions are related to the 
wave function and its derivatives at the origin x = 0, and the non-local contributions 
correspond to corrections to the wave function. The application of vacuum insertion to 
the present case of heavy- quarkonium decays is quite a strong assumption, because it 
removes already all the operators with a RCM derivative. This follows from the fact 
that, with this assumption, the c.m.s. frame of the hard annihilation process has to be 
the same as the one of the heavy quarkonium. Furthermore, this assumption leads to 
vanishing matrix elements for all colour operators at the small scale /i. And 
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finally, in such a picture the contributions of the purely gluonic piece (^) vanish (e.g. the 
term F^^) in 

If we nonetheless use vacuum insertion at a small scale /i we find 



(92) 



The first term is the one familiar from non-relativistic potential models: Identifying 

o!i^ = h'^m' = ^\m? (93) 

we recover the well-known Royen-Weisskopf formula (M = M(^Si) ~ 2m): 

T{'S,^a) = ^^ \Rm' , (94) 

which holds for the exclusive decay. However, it is this expression we obtain to leading 
order also for the inclusive decay; this indicates that the exclusive mode ip — >■ e~^e~ will 
saturate a large portion of the inclusive decay ip —>■ e+e~+ light hadrons. 

The non-logarithmic part of the Cg^'' correction in (^) has also been discussed 

by BBL 1^ and by Keung and Muzinich (KM) [l^. It is in fact the only correction that 
occurs in either approach. BBL propose to identify (the analogon to) c'q^ with the limit as 
r ^ of — A^i?(r) with appropriate regularization. KM calculate kinematical corrections 
to the leading-order amplitude. They denote by e the binding energy and by p the relative 
three-momentum of the heavy quark and antiquark. Then they evaluate the amplitude 
at = me rather than at = 0, as is the usual wave-function prescription for S-wave 
decays. Their corrections proportional to e/m are thus •p'^/im? corrections and can be 
identified with those arising from the ii^-operators in our approach. 

Let us emphasize that the appearance of the logarithm in ( |92D indicates the breakdown 
of the naive potential-model calculations also for S'-wave decays as was conjectured a long 
time ago |T^. Expanding 

Inr/ ^ ^ ' In — , (95) 



33 — 2n/ TT n 

we obtain a contribution to the decay width proportional to |i?"(0)p ln(m//i). This 
manifestly violates the potential-model ansatz that the infrared dynamics is given by the 
meson's non-relativistic wave function (and its derivatives), while the short-distance part 
is free of infrared singularities. 
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As a further example, let us study the hadronic decay of spin-triplet P-wave quarkonia 
n^Pj = xqj{i^P)- In the wave-function approach, the leading-order (in as{m)) decay is 
into two gluons for the J = 0, 2 states, r(^Po,2 99) = C'(a^), while the J = 1 state 
can first decay at 0(0^) into either three gluons or a quark-antiquark-pair plus a gluon. 



Moreover, it was found ||T9|, ^ |2T[ that the quark-antiquark-gluon cuts are singular in 



the limit of zero binding energy e 

,o , Hf 128 aAmY , m , , , ^,2 , ^ 

rePj^qqg) = ^ — ^^ln- \R[pm' . (96) 

The magnitude of the logarithm is usually estimated by identifying 1/e with the av- 
erage radius of the states. Again, the presence of an infrared sensitive logarithm 
signals the breakdown of the usual factorization assumption that is behind the wave- 
function approach. It has been argued that the decay into qqg has to be considered as 
being of the same perturbative order 0(q;^) as the two-gluon decay since the expresion 
as ln(m/e)i?'(0)^ has to be considered as a new nonperturbative parameter besides i?'(0)^ 
describing P-wave decays ||2^, |^ . 

The present formalism reproduces this result by generating the large logarithm through 
the renormalization group running. In fact, formally the decay into a pair of quarks is the 
dominant one, since it is of the same order in as as the two-gluon decay but logarithmically 
enhanced, see ( |102| ) below. To see how this comes about consider the qq decays of the 
spin-triplet P-wave states in more detail. At the matching scale m only colour-octet 
contributions are present (section 5) and hence the leading-order result is 

^ro(3x ^ qq^ light hadrons) = -^^a^(m)?2/(V'|Af = ^^al{ni)n f . (97) 
For the second-order contributions we find 

+c«''(Ef ,/i)(V'|£;fV), + c^''(Ef ,/x)(v^|4'V)m 

+C^^(Ef ,/i)(^|4'V)M + C'^\Af,ii){kf^ + + P(«)) . (98) 
The coefficients of the local operators are (at a scale fi < m): 
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C«9(Z)f ) = -l^al(m)nf Inr/ 

c«5m«) = — ^^2/ X J 

^ ^ ' 27m2 ^ -^33-2^/ ' 

mEf,,) ^ -J^^lim)2 (l - f ^m.) . (99) 

It is interesting to note that the renormahzation group flow induces operators {D^^ and 
D^^) which have not been present at the matching scale. We may express the local 
contributions to the second-order contribution in terms of the parameters bf'\ di^\ 
and ef^: 



126^^ fl — Inn 1 +8ef^ (l lnr7 ) 

^ V 33 - 2n/ 7 V 33 - 2nf ') 

128 ^(1) 1 gQ^(8) 1 _ ^(8) _ ^(8) _ ^(8) 

3 ^ 33-2n/ ' ^ 33-2n/ ' ^ 



.(100) 



We note that not only the leading contribution (^), but also the first subleading one (|100|) 
is the same for all three states ^Pq, ^Pi, and ^P2- Thus, if the channel involving a qq- 
pair really dominates the decays of these states into light hadrons due to the logarithmic 
enhancement, we have 

r(^P2 light hadrons) = T{^Pi light hadrons) = r(^Po ^ light hadrons) . (101) 

The second-order contribution ( |100| ) is given in terms of seven parameters, but this 
number reduces once vacuum insertion is assumed. Then only one of them is non-zero, 
namely the parameter di'\ and the rate takes the simple form 

327r 1 
r(\o m{9) light hadrons ) = -^al{m)nfd^l'^ ^^—^Inr] . (102) 

Furthermore, in a wave-function model the parameter d^l"* is given in terms of the deriva- 
tive of the wave function at the origin 

Sp{n) = A |Kp(0)|^ , (103) 

and if one expands ( |102| ) again in powers of as{m), see (p5|), we reproduce the result (pG]). 
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9 Comparison with previous approaches and conclu- 
sion 



Traditionally inclusive hadronic (and electromagnetic) decays of heavy quarkonia are 
calculated in the framework of the Bethe-Salpeter description of the bound state. As- 
suming an instantaneous potential and working in the extreme non-relativistic limit, the 
decay width of a heavy quark-antiquark bound state into light hadrons (l.h.) is written 
as 

r(n2^+^Lj ^ l.h.) = G{n) T{QQC^+^Lj) partons) . (104) 

This ansatz of separation into a short- distance part T describing the decay of a free (un- 
bound) QQ-pair to decay into partons (gg-pairs and gluons) and a long-distance one G{n) 
representing the non-perturbative bound-state formation, is motivated by the observation 
that the problem involves two widely separated scales, the QQ radius of the order of 1 fm 
and the heavy-quark Compton wavelength Aq ~ 1/m <^ ^qq- The non-perturbative part 
G{n) is expressed in terms of the non-relativistic (Schrodinger) wave function at zero rela- 
tive coordinate Rnsi^) (for S-wscve decays, and the derivative R'^p{0) for P-wave decays), 
which is calculated in a potential model or extracted from data. The factor G{n) thus 
depends explicitly on the binding energy e. The QQ decay, on the other hand, is governed 
by a scale of the order of ~ m, and can hence be expanded in a series of as{m). The 
decay rate is called factorizable into a long- and a short-distance contribution |jT9| if the 
short- distance factor T is calculable without encountering infrared divergencies. 

Let us elaborate on the assumptions of the conventional approach in some more detail. 
Consider the transition amplitude A4{J^,P; ki) of a quark-antiquark bound state of 
mass M (and four- momentum P = (M, 0)) into partons (light gg-pairs and gluons) of 
four-momenta in terms of which the width is given by r ~ \M\ dPS. (For illustration 
we suppress colour and polarization indices.) Introducing the Bethe-Salpeter bound- 
state wave function $(J^, P; g), dependent on the relative quark and antiquark momenta 
g {pq^q = |P ± g), the transition amplitude is written in this picture as 

M = I ^Tr <|.(P, g)0(P, g, (105) 

such that 0{P, g, k^) is the amplitude of a free (unbound) QQ-paii to decay into light 
hadrons. Upon a non-relativistic reduction 

<l>(P,g) = 27r5(g°) ^nUo) {ImSS, \J J-,) PssXP.q) , (106) 



35 



the amplitude becomes 



= E / S^^i^nUq) (ImSS, \JJ,)TiPssAP, iWP, h\ (107) 



m,Sz 



where pf = — = m? + em — g*^, e is the binding energy, and Pss^ is a projection 

operator. The Schrodinger wave function in coordinate space extends over distances of 
the order of the Bohr radius. Correspondingly, the Fourier transform ipnim{<l ) is non-zero 
only for \q\/m <^ 1, and hence Iql/m becomes a reasonable expansion parameter. For 
S'-wave decays, the leading term of this expansion yields the usual expression for the decay 
amplitude proportional to the non-relativistic wave function calculated at the origin in 
the relative coordinate space: 

M ^ -^RnsiO)A{q^ = 0, e) . (108) 



Here -R„i(0) = AiripniiO) {I = S,P,D, etc.) is the radial wave function at zero distance 
r = 0. For P-wave decays, the terms linear in q in the expansion of O, as well as terms 
linear in q coming from the small components of the relativistic wave function, must be 
retained yielding a final expression for the amplitude proportional to the derivative of 
the wave function for the / = P state at zero r, R'^p{0). Factorization is now said to hold 
if the limit e — > exists for the amplitude A. 

Within this potential-model approach, next-to-leading-order (NLO) perturbative QCD 
corrections have been calculated. Both the NLO corrections to ^5*0 (i.e. rjc) and ^Si 
(i.e. J/ip) decays indeed obey the factorized form (|104|) . However, there are three 



observations from which it becomes obvious that this picture is too simple: 

• Infrared sensitive logarithms ~ ln(m/e) appear in the calculation of P-wave de- 
cays, to be precise in the NLO corrections to ^Po,2 decays |20| and already in the 



leading-order expressions of ^Pi and ^Pi decays [0. That is, without keeping the 



binding energy non-zero, the perturbative part of the calculation would diverge. 
The factorization ( |104| ) thus breaks down. 

The description of S'-wave decays based on the strict non-relativistic limit is not in 
agreement with data. For example, a reasonable as determination from the ratio 
T(^Si l.h.)/r(^5'i — > ii) is possible only once a (rather large) adhoc relativistic 
correction factor is applied. Other failures are the photon spectrum of ^5*1 decays and 
photo- and hadroproduction of J/^/^ |]^. In fact, logarithmic infrared divergences. 
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predicted to arise also in relativistic corrections to S-wave decays a long time ago 
0, show up, cf. ( ]9^j95|) . Within the wave-function approach one can distinguish 



two sources of relativistic corrections, namely corrections to (i) the amplitude O and 
the (ii) the wave function ip, cf. ( |107|) . Relativistic corrections to the wave function 
are very difficult to access since at present no full analysis of the spin-independent 
relativistic corrections to the potential exists, cf. (^31). Kinematical corrections of 
type (i) have been discussed by Keung and Muzinich some time ago |]18| for S-wave 
decays and more recently applied also to J/ip photo- ||25| and hadroproduction 0. 
The prescription is to evaluate the amphtude A{q^,e), (|107|) , at = me rather 



than at q"^ = 0. In this way certain binding-energy corrections oc e/m = q"^ /m? are 
kept. 



• Sizeable, non-perturbative corrections have been predicted [^B] for annihilation de- 
cays of heavy quarkonia, at least for the charmonium system. Such higher-twist 
corrections could arise from non-zero condensates (in particular, the gluon conden- 
sate) and/or colour-octet intermediate states, e.g. via a non-perturbative colour-El 
transition of the J/ip into a (coloured) Xcj state ^Si —>■ ^Pj^^ + g, followed by a 
"hard" (/i ~ m) decay ^Pf'^ — > gg. 

In this paper we aim at a systematic, QCD-based treatment of inclusive annihilation 
decays of heavy quarkonia. In our approach, the factorization of long- and short- distance 
contributions is well defined. Furthermore, "genuine" relativistic corrections, i.e. the ones 
proportional to the relative velocity of the heavy quarks in the conventional language, 
and what is usually called non-perturbative corrections come out to have the same origin, 
namely the higher-order terms of the 1/m expansion. In fact, since we may shift certain 
contributions from the operators into the states in this expansion, there is no unique 
distinction between the two kinds of corrections. 

To achieve such a systematic description of quarkonia decays we formulate an effective 
field theory for heavy quarkonia, which follows rigorously from the QCD Lagrangian. As 
a first step we write down an operator product expansion for the inclusive annihilation 
decays of heavy quarkonia in order to separate long- and short-distance contributions. 
Herein the distance scale is set by the Compton wavelength of the heavy quark. This 
provides us with a systematic 1/m expansion of the short- distance contributions with 
mass- independent operators. As the next step we expand the heavy-quark fields Qlf^\x) 
and the Lagrangian. In the wave-function picture, the expansions of the states and of the 
kernel correspond to corrections to the amplitude A in (107), while time-ordered insertion 



of the higher-order Lagrangian correspond to corrections of the wave function. 
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In the set-up of an effective theory approach to heavy quarkonia we now observe a 
crucial difference compared to heavy-hght systems: The static hmit does not exist for a 
heavy quarkonium state. Divergent phases appear in the quark-antiquark sector of HQET 
once the velocities of the two heavy quarks differ by an amount of only Kqcd/^- These 
divergent phases can be and have to be absorbed into the heavy-quark states. Using 
reparametrization invariance, we have shown that it is necessary and sufficient to include 
the 1/m kinetic term into the leading-order Lagrangian in order to consider two heavy 
quarks moving with the same velocity. The dynamics defined by this static-plus-l/m- 
kinetic-energy Lagrangian already contains the divergent phases (and the binding of the 
two heavy quarks) as an infrared effect, and this is the physical reason why the limit 
V ^ v' exists. Since the redefined states of a heavy quarkonium (seemingly) do not have 
a static limit, it is, in contrast to heavy-hght systems, not possible to describe the full 
mass dependence of the relevant matrix elements in a 1/m expansion. Although the short- 
distance part may still be written as a 1/m expansion with mass-independent operators 
(once divergent imaginary parts have been shifted from the operators into the states by 
a suitable redefinition), the matrix elements of these operators with the redefined states 
become mass-dependent due to the mass dependence of the states. 

As a result of the mass dependence of the states, heavy-flavour symmetry does not hold 
anymore. However, judging from the spectra and also from the widths of charmonium 
compared to bottomonium, the flavour dependence of the matrix elements does not seem 
to be very strong; at least it appears to be much weaker than one would expect for a 
Coulombic system. 

On the other hand, spin symmetry is still preserved since the extra (kinetic) 1/m term 
in the leading-order Lagrangian is spin-symmetric. Hence the heavy quarkonia states 
described by our leading-order Lagrangian have to fall into degenerate spin symmetry 
quartets. Spin symmetry, furthermore, allows a restriction of the number of independent 
parameters describing matrix elements involving heavy quarkonia states. 

We have applied our approach to inclusive heavy-quarkonia decays into light quarks 
and leptons. The effective theory machinery allows us to calculate the logarithmic depen- 
dence on the heavy mass m by studying the renormalization of the operators mediating 
the decay. The matching of QCD to the effective theory is performed at the large scale m, 
where the coefficients of the operator are determined by comparing the QCD result with 
the effective theory. This determines the initial conditions of the operator coefficients. 
The renormalization group of the effective theory then allows the coefficients to be run 
down to some smaller scale fi. We now see how the factorization assumption of ( |104| ) is 
generalized in a proper QCD treatment: A given inclusive annihilation decay is, in gen- 
eral, the sum of various contributions, each of which is the product of a non-perturbative 
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contribution, parametrized as the matrix element of an operator renormalized at a scale 
/i, and a Wilson coefficient evaluated at the same scale. 

In addition to matrix elements of local operators also non-local contributions appear. 
From the conventional point of view these correspond to corrections to the wave functions 
of the states. Among these non-local terms one may identify spin-symmetry breaking cor- 
rections, which in the wave-function language generate differences between, for example, 
the ^Sq and the ^5*1 wave functions. 

The rate to leading order is given in terms of forward matrix elements of dimension- 
six operators, which at the one-loop level are scale-invariant, and hence no logarithmic 
enhancement through terms as ln(m//i) is present. Due to symmetries, the leading-order 
contribution is in general given in terms of only two parameters {a^l^ and af"^) for a given 
orbital angular momentum i (and a given n) of the heavy quarkonium state. 

The first subleading contribution is suppressed by two powers of the heavy mass and 
involves forward matrix elements of dimension-eight operators as well as non-local con- 
tributions corresponding to corrections to the states. The number of dimension-eight 
operators is quite large, although it is somewhat reduced by spin symmetry. In total, the 

(C) (C) (C) 

number of independent parameters is in general still large, namely 14 (6) , c} , d} , 
ef^\ k2^\ Cf^^ and F^*"^) for each heavy quarkonium angular momentum i. 

Intuitively one may expect that some of the parameters are smaller than others. If one 
assumes vacuum insertion, then all matrix elements may be interpreted as wave functions 
and its derivatives taken at the origin. In this way one is led to assume that all the colour 
combinations (g) T'* are suppressed compared to the 1 1 operators. Furthermore, the 
wave functions for P-wave states vanish at the origin, and this suggests that the matrix 
elements of the operators D!^^ dominate for the P-wave states. Let us emphasize, however, 
that our analysis is independent of these assumptions. Although our approach has many 
common features with the expansion advocated by BBL ||5[ it is more general since we 
give the expansion before additional, less rigorous assumptions have been applied: 

1. Within a set of operators of a given dimension we do not neglect those that are 
suppressed by powers of the relative QQ velocity v/cin the wave- function language. 
For instance, we do not neglect an operator involving a gluon field strength [iD, iD] 
in comparison to {iDy, although the two operators correspond to different powers of 
v/cin the language of 0. Such a procedure is adequate for Coulombic systems but 
may be a bad approximation for charmonium and bottomonium, which certainly do 
not behave Coulomb-like. Yet, up to order (A/m)^ no difference appears between the 
two approaches, since the matrix elements of antisymmetric products of covariant 
derivatives vanish due to spin symmetry. 
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2. We do not restrict the heavy quarkonium to the leading Fock state. Corresponding 
arguments are again based on counting of powers oiv/c: Using perturbation theory, 
each additional gluon associated with the (assumed) dominant QQ pair is ascribed 
an extra power of f/c via the identification v ~ asimv), valid for a colour- 
Coulomb potential. Although this estimate may be underlined by the multipole 
expansion, we do not know of any rigorous derivation. 

3. The non-local contributions from the time-ordered products with the Lagrangian in 
r2 ( |55D allows us to rigorously define the spin-symmetry breaking terms that appear 
at order 1/m?. Thus in principle we can calculate the coefficients that appear 
in relations as = R-q^[l + 0{v^/c^)], although in practise additional unknown 

(c) 

parameters enter, namely G] for each /. 

If we apply these additional assumptions, a large number of our non-perturbative param- 
eters may be dropped, and our approach yields the same result up to order (A/m)^ as the 
one of ref. p. 

In conclusion, the approach presented here is based on QCD and hence provides a 
model-independent basis for the description of heavy-quarkonium physics. It allows us to 
separate long and short distances, where the short- distance contribution may be evaluated 
perturbatively. The long-distance part is parametrized in terms of matrix elements, which 
involve up to order operators of dimension six and dimension eight. These matrix 

elements are beyond the effective theory approach and have to be taken either from data, 
estimated via non-perturbative methods such as sum rules, or they may eventually be 
calculated from lattice QCD. 
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APPENDIX 



A Spinorology 



Similar to the case of only upper components p^, we also have only four independent 



matrices once we project on upper and lower components. The projectors are 

P+ = \{l + i^) P+ = \{l-i^) (109) 

and we have the mapping of the sixteen Dirac matrices 

1 P+P- = 7m ^ P+1,P- 
7^.75 — ' P+'yf.-f^P- = P+l^P-Vf, 75 — > P+I5P- 
{-tWu P+i-t)a^uP- = P+K7. - Vu^^)P^ . (110) 

We chose the four matrices to be P+j^^P- (corresponding to the unit matrix) and P^'-fnP^ 
(corresponding to the three Pauli matrices). Note that v^P+'jf^P-. = 0, so these are indeed 
only three matrices. Hence, stricktly speaking, all Lorentz indices in ( ^ ) and (pT|-|5^) 
should be written as perpendicular indices only defined by 

< = ap {g'^^-v'^v^ . (Ill) 

Note that 

9'" = 3 . (112) 

One may obtain the projections for any Dirac matrix F in terms of these four matrices 
by the trace formula 

P+rP_ = i Tr {rP+75P-} P+75P- + ^ Tr {rP+7^P-} P+7^P- • (US) 

The projections with P_ and P+ interchanged are obtained by replacing v — > —v in the 
above equations. 

B Matching calculation to 0(l/m^) 

In this appendix we collect a few useful relations that appear frequently in the matching 
calculation. 



41 



B.l Expansions of bilinears 

Using ([25|), P^), and (|i3| ) we obtain, for the bilinears up to and including order 1/rn?: 
g(+)rQ(-) = u+)Th^-) 

+ _L 5^ (/,(+) [r, /,(-)) + |r, 2 ^ I j 



With the help of ( |113|) and the equations of motion iv-d'^ 

we find the following expansions of the bilinears up to and including order 1/m?: 

I m 

Q(+) 75 Q(-) =A^ + — ^ {-A, + 2A,2 + - A,r) 
lb m"^ 

Q{+)^«g{-) = 1 ^"^3 

2 m 

Q(+) 75 7° g(-) = + (A- + A?) 

2 m 

16 



q(+) (_,) ^-/3 g(-) =v''A§- A^ 



+-L (Af - Af - Af 
2 m ^ 



~'~16 m 
Here we have introduced 



^— {v^ - 2 + 2 - < + < - <) - (« /?)) • (115) 



Ai = 75 h^'^ 
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A'^ = /i(-) 
Al = 75 

A'f = ie^'^^PVf, h^+'^^^i /i(-) 

= Z 9^" i 9^ 75 /l(-) 

A^^o = z i 5^ 7^" /i(-) 

= i 9^" i d^p 7^ /i(-) 
Ai2 = 2e"^^^t;„^5^/iW7.^^ dI /i(-) 

= 2 e^^^P V, I 75 ^ 5^ /^^-^ 

Am = h'-^U D^P -f5i h^-^ 
A^^ = h^+U D^P -i^'^i h^-^ 

= i e'^PP'' I -f^i h^-'^ 

= z e'^PP^v^ h^+h 75 z /i^^^ . (116) 
B.2 Expansions of quatrilinears 

With the help of the bilinears (|116|) , we obtain 1/m expansions of quatrilinears Ai^{g,e}® 
A*. Owing to the fact that quatrilinears have to be parity-even, and that total derivatives 
may be neglected, we find that there are two dimension-six operators, no dimension- 
seven operator, and 17 dimension-eight operators. Three more dimension-eight operators 
(C3 -|- C4, + E2 — -E3) appear only if also the kernel is expanded. Including colour, 
there are twice as many operators 

= T''^T^O'-^\ (117) 
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The dimension- six operators are 



- A^^^ = AiAl 



A? 



AO. A* 



The dimension-eight operators that can be constructed from ( |116| ) are 



(1) 



D 



(1) 



(1 



-B. 



CI 



(1 



-a 



(1 



(1 



D, 
D 



(1 



(1 



El 
-E, 



Aa A* 

A3Al = 

AOtf^ A* 
^5 ^5af3 



= A,Al 

Aa A* 



Aa A* 
^W^2a 



1 



I e 



a/3fiu 



V,, A 



A* 

80/3 



2~ ~ ^ 
A^Al + A^Al 

le'^^^^'v^A^^pA*, 

Ai2A\ + AiAl^ = - 

Aa A* I A°' A* 



±AaP A* 
2^8 ^8ap 

A.At 



( A" A* -I- A"" A* ) 
A°'^ A* 



^8 ^5a/3 



Aa A* 



7a 



--le'^^'^'-v^AsafsAl + h.c 



AuAl + AiAl^ 
A°' A* -I- 4" 4* 

^16^2a + ^2^16a 

4° 4* _i_ 4° A* 

^15^2a + ^2^15a 

+ A^Al^ 

Aa A* I /!« 4* 
^18^2a + ^2^18a • 
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B.3 Two- component representation 

In the rest frame of the quarkonium we have v^^ = (1,0). Using the Dirac representation 
of the gamma-matrices let us introduce the two-component spinors ip and x via 



Then we can express the bihnears as follows 



A2 
A3 

A, 

A'i 

^6 

A'i 

Ax2 
^13 

^14 



-i d ■ [tl^^ ax) 

—'^^i D -ax 
— i ■j/'M (J X i D 



X 



^ - (id) 



-id {id - (V'^CTx)) 
i e'^^ i di aj i Dk x 
-I e'^'^i djijU dI X 

-^''(i d] X 



X 



(123) 
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-ip^ a [i D ] X 



Aie = -V'^ 



i D \i D -a 



An = i e'^^ Di aj i Dk x 



i D -a] i D 



X 
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Ak 
^18 



(124) 



The complex-conjugate expressions are obtained from ( |124| ) by the replacements V'^ ~^ 
X —>■ ip , i —>■ —i (using aj = ai). The quatrilinears are obtained using e^j^eiji = 2 6ki- The 
connection to the operators of ref. is then straightforward (c = 1 or 8) 



4 ^ 



D 



D 



4^^ 



-OXSo) 

-aC^i) 
a('Pi) 

3a('Po) 

a('Po) + a('Pi) + a('P2) 
a('Po)-a('Pi) + a('P2) 



(125) 



B.4 Matching calculation 

The short-distance coefficients appearing in (^9]) and (|55D at the scale fi = m are obtained 
from the annihilation part of the scattering amplitude Ai for QQ QQ computed in 
full QCD. The amplitude, calculated in QCD perturbation theory for on-shell quarks and 
antiquarks, is Taylor-expanded in the "small" components of the heavy quark momenta. 
The small component of a momentum is defined as the component perpendicular to the 
direction of the quarkonium momentum Mv. Denoting the momenta of the initial QQ 
pair by p and p for Q and Q, respectively, and the ones of the final QQ-pair by p' and p' 
for Q and Q, respectively we have: 



p = v\liv? — p\ + pi_ 



mv 



2m? 



{P + Txf + ... 



+ P + 7r 
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p = mv 
P' 
P' 



mv 



mv 



1 
1 

2m^ 



iP-nf + ... 
{P + r]f + ... 



+ P-n 
+ P + ri 

+ P-V 



(126) 



where vr^ = rj'^. Note that both the kernel and the spinors have to expanded to the desired 
order. The latter expansion follows from (EH). For example, up to order we have 



Q 



(+) 



1 + 



p2 + yr2 



1 + 



6m^ 



l + f 
2m 



P^ 



(127) 



The result of this procedure is then matched to the operators in ( ^91 ) and ([55| ) by identi- 
fying the momenta with the derivatives appearing in the operators 



(128) 



Then the coefficients C(Of ^) of can be read off from the dimension-eight contribution 
to the amplitude 



M 



d=8 



\2mJ 



E 



EC(i?f,m)i?r+5:C(Cr,m)C^ 



Ac) 



Ac) 



c=l,8 li=l 



+ Y.C{Dt\m)Dr + Y.C{Ef,m)E,, 



i=l 



Ac) 



(c) 



Ac) 



i=l 



i=l 



c=l 

+C(P:^-',m 



4m^ 

+ C(Ci'\m 
+ C(Cf,m 
+ C(Cf,m 
+ C(Cf ,m 
+ C(D 



75 <»75 



1 ? 



7o 



f ®f + C{Bf ,m) P"^ -i'^ 
(7r + r7)-P75®75 + C(C^'=\m) ®f +^ 
(tt + r/) • P7" ® 7„ + C(Cf\m) [f ® ^ + 
{-i) e"^^-^ t;^ P/3 K 7<5 ® 75 + V-y 75 ® 7^] 
{-i) e"^^"^ P/3 K 75 ® 7<5 + ^7 75 ® 75] 
TT ■ r/75 + C(P)^''\m)y (g)|f 



®?f ] 
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+ C{D^\m) {-i) e"^^^ vr^ r]^ [75 ® 75 + 75 ® 75] 

+ C(El^\m) [TT^ + r^^] 75®75 + {C(Ef\m)+C(Ef ,m)} ® ^ + ^ 

+ C(Ef [^r' + r/'] 7"® 7a} , 

where we suppressed the spinors and abbreviated the colour structure 

Xi = l®l , Xg = ® TV 

Note that and Eq give only field-strength terms, i.e. vanish for commuting 
matrices. The dimension-six part is simply 



(129) 

(130) 
Dirac 



M 



d=6 



c=l,8 1=1 



= A E {C(4'\"^)75®75 + C(4'\"^)7"®7a} • (131) 

^ c=l,8 

C Operators in terms of quarkonia quantum num- 
bers 

In order to calculate the hadronic matrix elements using vacuum insertion, it is convenient 
to re-express the operators B, C and D in terms of other operators in which for both of the 
heavy-quark bilinears the "orbital angular momentum" corresponding to the derivative 
is coupled to the total spin of the quarks to some total angular momentum J. If is 

either the RRM or the RCM derivative {i D ^ or id^, respectively), these couplings are 

i (/.(-)7^V^/i(+)) (/.(-)7'^V./i(+) 



'Pi 



7m 



-2 ^ 3p^ ^ ^hi-)^[^^^^]hM^ (/i(-)7[^V,]/iW) (132) 
where the bracket denotes the symmetric traceless combination of a tensor (cf. ( |111| , |112D ) 



X 



(133) 
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The first operator of tlie tliree sets B, C and D is already tlie ^Pi combination; tlie otlier 
tliree may be rewritten as linear combinations 







®'Po\ 


( 1/3 









^''Pi\ = 





1/2 




VP2 


®^P2) 


1-1/3 


1/2 


By inverting (|13^ 


Ep we 


have 













f'Po 
'Pi 
VP2 



3 p 



'Pi 



'P. 



Explicitly, for the D-operators we have, for example, 

^(c)^3p ^ _ 1 nW 



-D^2 
3 ^ 



2 



2 L - ■ - J 3 
Similarly we can define operators with appropriate S-wave quantum numbers 



El' 



2 

is 

2 

1 

4 



(c) 



(134) 



(135) 



(136) 



(137) 
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